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We study the electronic states of the anisotropic triangu- 
lar lattice Hubbard model at half filling, which is a simple ef- 
fective model for the organic superconducting k-BEDT-TTF 
compounds. We treat the effect of the Coulomb interaction 
by the fluctuation exchange (FLEX) method, and obtain the 
phase diagram of this model for various sets of parameters. 
It is shown that the d-wave superconductivity is realized in 
the wide region of the phase diagram, next to the antiferro- 
magnetic states. The obtained phase diagram explains the 
characters of the experimental results very well. 

PACS numbers: 74.70.Kn, 74.20.-z, 74.25.Dw, 75.50.Ee 



It is known that the layered organic compounds made 
of BEDT-TTF molecules exhibit rich variety of ground 
states, through the strong correlation effects between 
electrons. Among them, K-polytype of BEDT-TTF sys- 
tems has attracted much attention because many of them 
are superconducting materials. The superconducting 
transition temperature (T c ) reaches ~ 10K, whose ra- 
tio to the band width is as large as that of the high-T c 
cuprates. 

In fig. 1, we show the schematic pressure-temperature 
phase diagram of the k-(BEDT-TTF) compounds. For 
example, Cu[N(CN)2]Cl salt at ambient pressure is in the 
antiferromagnetic (AF) insulating phase as shown by the 
arrow in the figure. The Neel temperature (In) is 27K 
and the spin moment is greater than 0.4^b per magnetic 
unit, i.e., dimer of BEDT-TTF molecules. The antifer- 
romagnetism disappears with increasing pressure to give 
birth to the superconductivity. T c takes its maximum 
value, 13K, at 200bar, but this superconducting (SC) 
phase also suppressed by several kbar, where a param- 
agnetic metallic phase is stabilized. This means that the 
superconductivity coexists with the strong antiferromag- 
netic correlations. Thus, it is natural to consider that 
the antiferromagnetic fluctuation is the origin of super- 
conductivity. 

In the previous studies, the organic compound has 
been studied theoretically through the local density ap- 
proximation and the Hartree-Fock approximation. 
The latter study could explain the pressure dependence 
of the AF boundary in fig. [j]. The structures of these 
compounds are very complex. Nonetheless, the electronic 
structure near the Fermi level is very simple, because it is 
composed of well separated HOMO levels. The simplest 
effective model is a triangular lattice Hubbard model 



with anisotropic hoppings at half filling, as shown below. 
The effective on-site Coulomb repulsion (U) is estimated 
to be comparable with the band width (W) in this model. 

The aim of this letter is to discuss the mechanism of the 
superconductivity, and to understand the entire phase 
diagram shown in fig. 1. For this purpose, we study the 
anisotropic triangular lattice Hubbard model by using 
the fluctuation exchange (FLEX) method. It is one of 
the methods that have much advantage for the systems 
with large spin fluctuations and that can handle the SC 
and AF states on the same footing. In this letter, we 
find the d-wave SC phase next to the AF phase. The 
obtained T c is consistent with experiments. This is the 
first systematic study of the organic SC compounds by 
using the FLEX method. 

The FLEX method is a kind of self-consistent pertur- 
bation theory with respect to U, which gives reliable re- 
sults although it is an approximation. For example, the 
imaginary time Green function obtained by the FLEX 
method agrees well with that by the QMC for the square 
lattice Hubbard model with moderate U. B This method 
has been applied to the study of high-T c , and various 
non-Fermi liquid behaviors are reproduced well. |^ [ll]] 
Recently it has also been applied to the superconducting 
ladder compound, Sri4_ a: Ca x Cu2404i. [O] 

Dimerization of k-(BEDT-TTF) salts is so large that 
we can take a dimer of BEDT-TTF molecules as a unit. 
So, as a simple effective model, we consider the two- 
dimensional single band Hubbard model with the lattice 
structure shown in fig.|[ where to, t\, t<x are hopping pa- 
rameters. We also take into account the on-site Coulomb 
repulsion (U), which is an energy loss when two elec- 
trons are on the same BEDT-TTF dimer in the original 
k- (BEDT-TTF) lattice. The hopping parameters depend 
on materials, but both the band calculations and experi- 
ments indicate that t\ ~ to and |ta| -C |to|, [213- 1^] 
In this study, we put to = — lj ^2 = —0.05, and change 
ii(j0) as a parameter. In this sense, this model is an 
anisotropic triangular lattice. For |ti| < |io|, W is al- 
ways 8 1 to | . A similar model has been studied recently, 
but the essence of the physics of k-(BEDT-TTF) salts is 
contained in the above simple model. |l(J 

Here, we explain the FLEX method. The Dyson equa- 
tion is written as 



{G(k, e,,)}- 1 = {G°(k, e n )} 1 - S(k, e n ), 



(1) 



where G°(k,e) is the unperturbed Green function. The 
self-energy is given by 
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E(k,e„) =T5^G(k-q I c B -« I ) -U 2 
1,1 

X Q X (-> (q, cut) + i X W (q, un) - X°(q, , (2) 

X ( ± )(q^)=X -{i±C/x (q,^)}^, (3) 
X°(q,w,) = -T^G(q + k,^+e„)G(k, e „), (4) 

k,n 

where e n = (2n + l)7rT and ui = 2/7rT, respectively. We 
solve the equations self-consistently, choosing the 

chemical potential /j, so as to keep the system at half- 
filling. 

To determine Tjy, we calculate the Stoner factor with- 
out vertex corrections, as, given by 

a s = max{ [7- x (k, W = 0) }, (5) 

The antiferromagnetic critical points are determined by 
the Stoner criterion, as = 1. 

We also determine T c by solving the linearized Eliash- 
berg equation with respect to the singlet-pairing order 
parameter, tf>(— k, e n ) = +4>(k, e n ), 

A • 0(k, e n ) = -T^2 ^( k - q> e « - e ™) 

q , m 

xG(q, e m )G(-q, -e m ) • </>(q, e m ), (6) 
K(k,«0 = -tfV^Ck,^) - i(7 2 x (+) (k,^) + 17, (7) 

where T c is given by the condition that A = 1. 

The theories by Mermin, Wagner and Hohenberg pro- 
hibits finite 7n and T c in two dimension. |l7],[l8| It is well 
known that as given by eq.(||) satisfies this condition, be- 
cause FLEX treats the spin fluctuations self-consistently. 

So we determine Tn by the condition as = Qn, 
where we set aN as (1 — on) -1 ~ 0(100). The AF state 
will occur through the weak coupling between layers. On 
the other hand, A = 1 is fulfilled at finite T c by using 
eq.(|^). The obtained T c is, however, reliable in many 
cases. 

First, let us see the phase diagram as a function of 
1*1 1 (= -*i), U and T, for T > 0.02. A number of sec- 
tions of the phase diagram are plotted in fig.|[ 16x16 
/c-points and 256 Matsubara frequencies (abbreviated as 
16x16x128 mesh) are used to draw this phase diagram. 
For 32x32x1024 or 64x64x512 mesh the Stoner factor 
becomes larger, but that effect can be compensated by 
taking slightly larger aN- In fig. [| we set qn = 0.99. As 
for the superconductivity, the size effect of T c is rather 
smaller than that of Tn, because A rapidly grows near 
the critical temperature. 

Fig.||(a) shows the critical temperatures, Tn (small cir- 
cles) and T c (large circles), on the planes of [7 vs T at 
several values of |*i|(= — *i). Only AF phase appears 
at |ti |=0.3, SC occurs for |ti| = 0.4 ~ 0.7 on the left 
hand side of the AF phases. As |ii| increases, we see 



the followings: The onset value of U for the AF region 
shifts larger, at the same time, the peak of Tn decreases, 
and the SC region are extended. The maximum value 
of T c at each section increases for |ii| < 0.6, and then 
decreases for \t\ \ > 0.6. 

Fig.||(b) shows T c and Tn on the planes of |fi| vs T, 
at [7=3.0, 6.0 and 9.0 for T > 0.02. At U = 3.0 only 
the AF phase appears. On the other hand, the SC state 
occurs next to the AF phases at [7=6.0 and 9.0. As U 
increases, both the SC and AF boundaries move toward 
larger \ti\ region. At [7=9.0, which is the largest [7 in 
this figure, Tn decreases as |ii| becomes larger and the 
SC state takes place instead of the AF state for \ti | >0.7. 
This SC phase , however, disappears for |ii| >0.75. At 
the \ti |=0.8, the AF occurs at larger [7, but no SC state 
is stabilized at T > 0.02. No SC and AF phases appear at 
\ti |=0.9 and 1.0 for [7 < 40. We will discuss the reason 
for this in the latter section. The value of A increases 
monotonically not only as T is decreased, but also as 
decreasing \ti \ as shown in fig.|| 

Next in fig.|], we plot the SC order parameter, (p(k, ui = 
0), along the line with the arrows on the Fermi surface. 
The parameters used in fig.|] are [7=4.6, ti=— 0.5 and 
T = 0.024, where A = 1.02. One can clearly see that 
there are nodes in the (7r,7r) and (— 7r,7r) directions and 
the superconductivity has c? :c 2_ 1) 2-like symmetry The 
maximum values of 4>(kp, w — 0) is slightly shifted from 
(0,7r) direction to (— ir,ir) direction. This is also the 
case for all the other superconducting phases. We have 
checked that the odd-parity superconductivity is not re- 
alized in the present model. 

To understand the overall aspects of the phase dia- 
gram, we discuss it on the basis of the weak coupling 
theory. In this scheme, the effective interaction between 
a singlet Cooper pair is given by 

Aw = - max(y(g - q', O^WgOWd^HFS, (8) 

<t> 

where 4>{q) is the energy-independent SC order parame- 
ter, and V(q, uj) is given by eq.(0). The obtained 4>(q) has 
the d-wave-like symmetry in this model as shown above. 
Introducing the cutoff energy lu c , the T c within the weak 
coupling theory is given by T c = 1.14w c exp(— 1/A W ). 

In fig.[| we show the obtained x 1 ' {l^ = 0) by the 
FLEX method for |ii| = 0.3 - 0.7. At each |ii|, we 
choose [7 so as to satisfy as — 0.99. Then, V(q) w 
W 2 /2-x i ~ ) (q,0) and max q {V(q)} « 317/2(1 - a s ). For 
each X (<?! 0) shows a sharp peak around (n, it), 
which favors the d-wave SC state similarly to the high-T c 
cuprates. We note that x^ - "* (<7, 0) has an incommensurate 
structure when as is rather small. 

As \ti | increases from 0.3 to 0.7, the peak of KSjO) 
around (71", ir) becomes broader, because the nesting con- 
dition becomes worse. This situation will make T c lower 
within the weak coupling scheme. On the other hand, as 
|ti I increases, the value of ma,x q {V(q)} becomes larger, 
which will make T c higher. The optimum |£i| for SC, 
which is |ii| = 0.6 in our calculation, may be determined 
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by the balance of these two opposite effects. We have 
already shown that d-wave SC occurs at T > 0.02 in 
the case of \ti\ = 0.4 - 0.7. We also find T c w 0.01 at 
1^1 = 0.8. 

At \ti\ = 0.3, AF phase completely prevails over the SC 
phase because max q {V(q)} is rather small at as = 0.99. 
On the other hand, when \ti | > 0.9, as = 0.99 is not sat- 
isfied for U < 40. So, ma,x q {V(q)} is also small. Besides, 
the structure of V(q) changes drastically as shown in the 
inset of fig.[|.: The peak position of [Qj 0) shifts from 
(it, 7r) towards (27r/3, 2tt/3), so the AF correlation is vio- 
lated. If we map our model to the corresponding Hciscn- 
berg spin system on the anisotropic triangular lattice, 
the ground state classically changes from AF order to 
the 120 degree structure at Ji/Jo = (ii/io) 2 = 2/3, i.e., 
ti/to = 0.82. The momentum (2tt/3,2tt/3) corresponds 
to the 120 degree structure. As a result, for |ti| > 0.9, 
any singlet SC state is not realized up to U < 40. 

The Fermi surfaces at ti = —0.5 and T — 0.024 are 
drawn in fig.||. The change of the Fermi surface is at- 
tributed to the renormalization of \t\\ by many body ef- 
fect. The sum rule, which must be satisfied in principle, 
but not necessarily fulfilled in actual numerical calcula- 
tions, is checked to be satisfied within 0.5%. 

Let us discuss the connection between the real materi- 
als and the obtained phase diagram. The T c of k-(BEDT- 
TTF) is as high as 10K. If one takes the noninteracting 
band width (W) to be 5000K, then 0.03|t |=18K, which 
corresponds well with the experiments. H 

Next we consider the effect of pressure. First we return 
to the original k-BEDT-TTF crystal. It is natural to as- 
sume that the effect of pressure is most sensitive to the 
loosely-packed parts, i.e., to the interdimer ones. The in- 
terdimer transfer integrals determine the band structure. 
So t's or W in the present model is sensitive to the exter- 
nal pressure, while the interdimer transfer integrals and 
the on-site Coulomb energy on the BEDT-TTF molecule, 
or U in the effective model, is insensitive to the applied 
pressure. The ratio of ti/to will change, but both to and 
t\ are estimated from the interdmer transfer integrals, 
thus the ratio does not change so large as U /W. Because 
W is fixed in fig.||, one can regard the pressure effect as 
decreasing U. 

Some band calculations suggest that \t\\ ~ \tp\, where 
no SC and AF phases appear in fig|| however. |D| One 
possibility for this discrepancy is that the band calcula- 
tions estimate the value of tx/to larger. The other is 
that the current model is too simple to reflect multi- 
band effects properly. Such a effect is probably taken 
into account by the renormalization of t\/to of this ef- 
fective model, because the AF phase will emerges at 
Ji/Jo ~ 2/3, or fi/to < 0.82 in the present model as 
noted above. If one choose such a ti/to, then one can 
find not only the AF state, but also the SC state next to 
the AF phase. Moreover the SC phase appears in higher 
pressure region. In this sense, the phase diagram in fig.^ 
agrees well with the experimental phase diagram. 

Furthermore, the experimental AF state is commensu- 



rate, [Q which is also consistent with the FLEX results 
that the peak of x(l>fy is always at (tt,tt) in the AF 
phase, due to the self-consistency. Although the system 
is half-filling, the region where SC state appears in fig.|| 
will not be replaced by the nonmagnetic Mott insulat- 
ing one, because U is smaller than W in that region. 
We finally make a comment for chemistries that T c takes 
the maximum value, ~ 0.029, at t\ ~ 0.6. T c possibly 
becomes higher than that of the currently known ma- 
terials if one can compose such compounds, though the 
correspondency between the present simple model and 
the original complex structure is not straightforward. 

In conclusion, we have studied the anisotropic trian- 
gular lattice Hubbard models at half-filling by using the 
FLEX method. We obtained the phase diagram, where 
the ci-wave superconducting phase is realized next to the 
antiferromagnetic one. The resultant phase diagram is 
consistent with the experimental phase diagram in k- 
(BEDT-TTF) systems. These results suggest that the 
electron-electron correlation is important in /{-(BEDT- 
TTF) compounds and that the spin fluctuation mecha- 
nism is the possible origin of the superconductivity. 
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FIG. 1. Schematic phase diagram of k-(BEDT-TTF) com- 
pounds. 

FIG. 2. A model and the Fermi surface when ti = —0.5 
at T=0.024. The solid line shows the one at [7=4.6 and the 
dotted line, (7=0 . Arrows on the solid line are explained 
in fig.4. 64x64 fc-points and 512 Matsubara frequencies are 
used. 

FIG. 3. Calculated phase diagram on the planes of U vs T 
at several values of |ti| (a), and on the planes of |ti| and T at 
a few U (b), respectively. The antiferromagnetism is realized 
in the temperature region below the small circles, and the 
superconductivity, below the large circles. Lines are guides 
for the eye. 

FIG. 4. The momentum dependence of 4>{kF,u> = 0) along 
the line with the arrows on the Fermi surface in fig. 2 



FIG. 5. X (_) («> w ) from <7=(0,0) to (n,n) at ti = -0.3, 
T = 0.03 and U = 2.85 (solid line), ti = -0.5, T = 0.03 and 
U = 5.00 (broken line) and ti = -0.7, T = 0.03 and U = 9.60 
(dotted line). max fe Ux°(k,u> = 0) is around 0.992 at fe=(7r,7r) 
for all the cases. The inset shows X^ _ H?i a; ) fr° m 9— (0,0) to 
(n,n) at ti = -1.0, t 2 = 0, T = 0.025 and U = 15 (solid line) 
and at ti = -1.0, t 2 = 0.05, T = 0.02 and U = 15 (broken 
line), max, U^{q,uj = 0) at ti = —1.0 is about 0.953 at 
g=(27r/3,27r/3). The 64x64x512 meshes are used. 
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(-0.8671,0) (0,0.86ti) (0.86tc,0) 

Fig. 4 
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